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Abstract 

Context. Reconstructing the evolution of post-common-envelope binaries (PCEBs) consisting of a white dwarf and a main-sequence 
star can constrain current prescriptions of common-envelope (CE) evolution. This potential could so far not be fully exploited due 
to the small number of known systems and the inhomogeneity of the sample. Recent extensive follow-up observations of white 
dwarf/main-sequence binaries identified by the Sloan Digital Sky Survey (SDSS) paved the way for a better understanding of CE 
evolution. 

Aims. Analyzing the new sample of PCEBs we derive constraints on one of the most important parameters in the field of close compact 
binary formation, i.e. the CE efficiency a. 

Methods. After reconstructing the post-CE evolution and based on fits to stellar evolution calculations as well as a parametrized energy 
equation for CE evolution, we determine the possible evolutionary histories of the observed PCEBs. In contrast to most previous 
attempts we incorporate realistic approximations of the binding energy parameter A. Each reconstructed CE history corresponds to 
a certain value of the mass of the white dwarf progenitor and - more importantly - the CE efficiency a. We also reconstruct CE 
evolution replacing the classical energy equation with a scaled angular momentum equation and compare the results obtained with 
both algorithms. 

Results. We find that all PCEBs in our sample can be reconstructed with the energy equation if the internal energy of the envelope 
is included. Although most individual systems have solutions for a broad range of values for a, only for a = 0.2 - 0.3 do we find 
simultaneous solutions for all PCEBs in our sample. If we adjust a to this range of values, the values of the angular momentum 
parameter y cluster in a small range of values. In contrast if we fix y to a small range of values that allows us to reconstruct all our 
systems, the possible ranges of values for a remains broad for individual systems. 

Conclusions. The classical parametrized energy equation seems to be an appropriate prescription of CE evolution and turns out to 
constrain the outcome of the CE evolution much more than the alternative angular momentum equation. If there is a universal value 
of the CE efficiency, it should be in the range of a = 0.2 - 0.3. We do not find any indications for a dependence of a on the mass of 
the secondary star or the final orbital period. 

Key words, binaries: close - stars:evolution - white dwarfs 

1. Introduction to a CE configuration: the core of the giant (the future white 

dwarf) and the initially less massive (hereafter the secondary) 

Virtually all compact binaries ranging from low-mass X-ray bi- ^^^^ ^^^^^^ towards their center of mass while accelerating and 

naries to double degenerates or pre-cataclysmic variables (pre- fl^^jiy expelling the gaseous envelope around them. 
5—1 ' CVs) form through common-envelope (CE) evolution. A CE 

5^ ■ phase is believed to be initiated by dynamically unstable mass Although the basic ideas of CE evolution have been outlined 

ti-ansfer from the ev olving more massive star to the less massive already 30 years ago, it is still the least understood phase of close 

main-sequence star (PaczyriskilgTl |Webbinj| 1984t |Hjellming| compact binary evolution. Theoretical simulations have shown 

[1989) . This situation occurs especially if the evolving more mas- that the CE phase is probably very short, ;10^ yrs, that the spi- 

sive star fills its Roche-lobe when it has a deep convective enve- raling in starts rapidly after the onset of the CE phase, and that 

lope (usually on the giant or asymptotic giant branch). Then the the expected shape of post-CE planetary nebula is bipolar-. For 

radius of the mass donor may increase (or stay constant) as a re- recent theoretical models of the CE phase see Taam & Rickei 

sponse to the mass transfer, while its Roche-radius is decreasing. (120061) and references therein. Despite the central importance of 

The resulting runaway mass transfer drives the mass gainer out CE evolution for a range of astrophysical contexts, hydrodynam- 

of thermal equilibrium because it accretes on a time scale faster ical simulations that properly follow the entire CE evolution ai-e 

than its thermal time scale. Consequently, the lower-mass star cuiTently not available. Instead, simple equations relating the to- 

also expands until it also fills its Roche-lobe, which then leads tal energy or angular momentum of the binary before and af- 

ter the CE phase are generally used to predict the outcome of 

Send offprint requests to: M. Zorotovic CE evolution. These equations are mainly used with the struc- 
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tural binding energy parameter (A), the CE efficiency (a), or the 
angular momentum parameter (y), which are all treated as di- 
mensionless parameters. The numerical values of these crucial 
parameters have so far not been constrained, neither observa- 
tionally nor theoretically. 



2. The sample 

Our sample of PCEBs consists on 35 new systems identified with 
the Sloan Digital Sky Survey (SDSS) and 25 previously known 
systems. To obtain a homogenous sample of systems we ex- 
cluded several PCEBs that appear in previously published lists. 



iNelemans et aP (l2000h and lNelemans & Tout! (I2005L herafter 
NT05) developed an algorithm to reconstruct the CE phase for 
observed white dwarf (WD) binaries. They derive the possible 
masses and radii of the progenitors of the W Ds in these binarie s 
from fits to detailed stellar evolution models dHurlev et al.l200Ci|) . 
This information can then be used to recons truct the mass- 
transfe r phase in which the WD was formed. INelemans et alJ 
( 1200 0*) used this method to reconstruct the CE phase of dou- 
ble WDs and find that reconstructing the first CE phase of vir- 
tually all double WDs requires a physically unrealistic high (or 
even negative) efficiency. Later NT05 extended their analysis to 
PCEBs and found no solution for two long orbital period PCEBs 
(AYCet, Porb = 56.80 d; Sanders 1040, Porb = 42.83 d). This led 
the authors to the conclusion that the energy equation fails in ex- 
plaining CE evolution. They proposed to use angular momentum 
conservation instead because they find the predictions of this re- 
lation to agree with the properties of observed binary samples. 
As mentioned above, the proposed angular momentum equation 
is scaled with the y parameter, and NT05 show that the values 
required to reconstruct the CE evolution of close WD binaries 
cluster in the range of y ~ 1.5 - 1.75, which has been inter- 
preted as a strong argume nt in favour of the y- algorithm. Later 
Ivan der Sluv s et af] (l2006t) extended the study of'Nelema ns et alj 
(2000), including more double WDs and calculating the binding 
energy of the hydrogen envelope instead of assuming a constant 
value for A. Exploring several options and combinations for the 
two episodes of mass transfer they find that indeed the evolution- 
ary history of the observed double WDs cannot be reconstructed 
by two CE ph ases described b y energy conservation. However, 
more recentlv lWebbiiig (l2008h " showed that the evolution of the 
observed double WDs can be understood within the energy pre- 
scription if quasi-conservative mass transfer for the first phase of 
mass transfer, and mass loss prior to the second phase of mass 
transfer (the CE phase) is assumed. In addition, according to 
IWebbinlJ (|2008|) the two problematic long orbital period systems 
in NT05 are probably post-Algol systems, i.e. also the product 
of qua si-conservative mass transfer, and not PCEBs. Webbink 
(l2008h convincingly demonstrates that the internal energy of the 
env elope has t o be taken i nto account, as suggested earlier by 
e.g. lHan et all ([1991 l2002h in the context of extreme horizontal 
branch stars. 



In any case, it is important to keep in mind that all the studies 
of CE evolution mentioned above are based on the analysis of 
small and not necessarily representative samples of PCEBs. 
We are caracterizing the first large and well defined sample 
of PCEBs (Giinsicke et al. 2010, MNRAS in prep.) based on 
intensive follow-up observations of white dwarf/main-sequence 
(WDMS) binary stars identified by the Sloan Digital Sky 
Survev ([Rebassa-Mansergas et al. 2007; Schreiber et al. 2008; 
Rebassa-M anserg as et alJ i2008i; Nebot Go mez -Moran et aL 



20091; iPvrzas et al.1 120091; iRebassa-M ansergas et al.1 l2010l) . In 



2.1. SDSS systems 

The theoretical research presented here has become possible due 
to considerabl e observational efforts in the last decade. First of 
all, the SDSS dAdelman-McCarthv et al] 12008'; 'Abazaiian et aL 
|2009) prove d to efficiently identify WDMS st ars. Schreiber etaL 
(l2007l) and iRebassa-Mansergas et al.l (l2010i) presented com- 
plementary samples of ~ 300 and ~ 1600 WDMS binaries 
from the SDSS. We initiated an extensive follow-up program 
of these stars to identify and characterize a large sample of 
WDMS binaries that underwent CE evolution. The first observa- 



tional res ults have been p re sente d by Rebassa-Manserg as et al, 
(IQQI); 'Schreiber et al.' ('2008); Rebassa-Mansergas et aL 
(20 08.);. Nebot Gomez-Moran et al. (2009); Pyrzas et al. (200C 
Sc hwope et all (120091) and IRebassa-M ansergas et al. (2010). At 
the time of writing (March, 2010), we have measured orbital 
periods for 53 SDSS PCEBs. From this sample, we excluded 
PCEBs with DC/DB primary stars because reliable estimates 
of the WD masses are not available for these systems. We also 
excluded systems with WD temperatures below 12000 K if the 
parameters w ere determined by spect ral fitting methods. As 
mentioned by iDeGennaro et al.l (l2008l) . it seems that spectral 
fitting methods probably lead to systematically overestimating 
the WD masses of these systems. We kept eclipsing systems 
with WD temperatures below 12000 K (e.g. SDSS 1548 + 4057) 
because independent tests for the WD mass are available for 
these systems. In summary, we have reliable measurements of 
both stellar masses and the WD temperature for 35 of the 53 
SDSS PCEBs with known orbital periods. These 35 PCEBs 
certainly form the most homogeneous sample of close com- 
pact binaries currently available, and the observational biases 
affecting this sample are expected to be small, as discussed in 
detail in Gansicke et al. (2010, MNRAS in prep.). The new 35 
systems with reliable orbital parameters from SDSS are listed in 
TablelAH 



2.2. Non-SDSS PCEBs 



this paper we reconstruct the evolution of the new, large and 
more homogeneous sample of 60 PCEBs with the aim to derive 
improved constraints on current theories of CE evolution in 
general and the CE efficiency in particular. 



Based on Table Al in NT05, Table 1 in ISchreiber & Gansickd 
(2003^ and with some additional recent identifications from 
Burleigh et al. (2006');'Tappert et al.1 (l2007l) ; lDrake et all (l2009t) ; 
Tappert et al.. (20Q9) we compiled a list of PCEBs that were 
not identified with our SDSS PCEB survey. In order to ob- 
tain a homogeneous sample that contains only WDMS systems 
we excluded all systems with hot sub-dwarf primaries (KV Vel, 
MT Ser, NY Vir, HS 0705-h6700, PN A66 65, V477 Lyr, TW Crv, 
UU Sge, AA Dor, HW Vir). We also excluded all systems where 
either the orbital period, one of the stellar masses or the 
WD temperature was not measured properly (HS 1136-1-6646, 
G1781A, HD33959C, G203-047ab, V651Mon, BPM71214). 
For four additional systems observational results pointing to- 
wards a peculiar evolutionary history appeared in the litera- 
ture: Sanders 1040 and AY Cet have extremely low WD masses 
and are almost certainly post-Algol binaries inst ead of PCEBs 
as me ntioned by Webbink (2008). According to lO'Brien et al.l 
( 1200 ll) . the primary in V471Tau is probably the result of a 
merger (a blue straggler), so the evolution of this star cannot 
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be approximated by single-star evol ution. Another system we 
excluded is EC 13471-1258, because O'Donoghue et alJ (l2003h 
show that it is probably a hibernating CV instead of a PCEB. 
Our final set of 25 non-SDSS PCEBs is listed in Table [Q 



3. Post-CE evolution 



In this section we follow ^S chreiber & Gansickd (l2003h and re- 
construct the post-CE evolution of the PCEBs in our sample. We 
assumed two different prescriptions of disrupted magnetic brak- 
ing. The reason for the choice of disrupted magnetic braking is 
the convincing support of this hypothesis from observations of 
CVs: (1) Disrupted magnetic braking explains the famous orbital 
period gap, i.e. the significant deficit of CVs in the orbital period 
range of 2 - 3 h.; (2) the current mass-transfer rates derived from 
observations of CVs above the gap are significantly higher than 
those of CVs below the gap; (3) the mean accretion rates de- 
rived from accretion-induced compressional heating are system- 
atical l y higher above than below the gap (Townslev & Bildste^ 
120031: [Townslev & Gan sicke 2009); (4) the donor stars in CVs 
above the gap seem to be slightly expanded compared to main- 
sequence stars, which is consistent with the dono r stars being 
driven out of thermal equilibrium ('Knigge"2006'); and (5) we 
find the fraction of PCEBs among WDMS binaries to be sig- 
nificantly decreasing towards higher masses at the fully convec- 
tive boundary tSchreiber et al]|2010h which has been pred icted 
by disrupted magnetic braking ( Politano & Weilerl l2006l) . We 
here consider two forms of disrupted magnetic braking, i.e. clas- 
sical d isrupted magnetic braking according to Rappaport et al. 
(Il983h and a more recently developed prescription taking into 
account the expected decrease of magnetic braking when the 
size of the convecti ve envelope of the secon dary star de c reases 
(iHurlev et al.ll2002h . We furthermore follow 'Davis et all (l2008h 
and normalize the latter prescriptions to obtain agreement with 
the mass-accretion rates derived from observations of CVs above 
the orbital period gap. 

The next key ingredient for analyzing the post-CE evo- 
lution is to derive the age of the PCEBs by interpolat- 
ing cooling tracks of WDs. We used the cooling tracks by 
Althaus & Benvenutol Sl99% for He WDs (Mwd;0.5Mo) and 
WoodI (119951) for CO WDs (Mwd;0.5Mo). We then calculated 



the orbital periods the PCEBs had at the end of the CE phase 
(f ce)- The required equations fo r classical magnetic braking 
are giv en in lSchreiber & Gansickd (|2003in For the Hurley et al.l 
(l2002h prescriptio n of disrupted magnetic braking normalized by 
iDavis et all (|2008|) we obtain 



PcE = 



(3Ct, 



1 1 1 



(1) 



with fcool being the cooling age of the WD. The masses and the 
radius of the secondary are in solar units, the period in years 
and C = 3.692 x lO""". The mass of the secondary's convective 
envelope M2,e is given by 



M2,e = 0.35 



1.25 -M2 
09 



(2) 



for 0.35 < M2 < 1.25 (see'Hurle v et"al]l2000l) . 

In Tables lATI and IA.2I we list the stellar and binary param- 
eters of the PCEBs in our sample as well as their cooling age 



' W e found a typographic error in Eq. (11) in lSchreiber & Gansickd 
j2003h : 9n should be replaced by In. 



(fcool) and the orbital period the PCEB had at the end of the CE 
phase (Pce)- The corresponding orbital period distributions are 
shown in Fig.[T] As most of the observed PCEBs are relatively 
young and most of our PCEBs have low-mass secondary stars 
for which gravitational radiation is assumed to be the only sink 
of angular momentum, the reconstructed zero age post-CE dis- 
tribution of orbital periods is not dramatically different from the 
observed distribution. In addition, the distributions of the sys- 
tematically identified SDSS PCEBs (black histogram in Fig.[Tl) 
do not differ significantly from the distribution of previously 
kno wn PCEBs that have been id entified through various chan- 
nels jSchreiber & Gansickd l2003). In the following sections we 
use the zero-age PCEB parameter reconstructed with the dis- 
rupted magnetic braking prescription as giv en by iHurlev et al.l 
( l2002h and normalized bv lDavis et alj (12008 ^. After reconstruct- 
ing the post-CE evolution, we can now concentrate on discussing 
implications for theories of CE evolution that can be drawn from 
our sample. 



4. CE equations 

It is generally assumed that the outcome of the CE phase can 
be approximated by equating the binding energy of the envelope 
and the change in orbital energy, and by scaling this equation 
with an efficiency a, i.e. 



(3) 



where Eg, is the gravitational (or binding) energy and AiSorb = 
£^orb,i - ^orb.f IS the total change in orbital energy during the CE 
phase. A variety of slightly different expressions for iSorb.i, £^oib,f, 
and Egr appeared in the literature and we briefly review them 
here. 

The final orbital energy Eoih,t is always calculated as the or- 
bital energy between the core of the primary (Mi c) and the sec- 
ondary (M2) at the final separation (of ) 



■^orb,f 



1 GMi,eM2 

2 flf 



(4) 



In contrast, different descriptions exist for the gra vitational 
energ y and the initial orbital energy. Some a uthors (e.g. lWebbinkI 
IT984t Ide Kooll \l99dt. iPodsiadlo wski et all I200I calculate the 
gravitational energy as being between the envelope mass (Mi ^) 
and the mass of the primary (Mi) 



p = 



GM,M 



AR 



(5) 



where A depends on the structure of the primary star, and the 
initial orbital energy as the orbital energy between the primary 
and the secondary at the initial separation (ai) 



_ 1 GM1M2 

^orb.i ~ 

2 a\ 



(6) 



As in 'Kiel & Hurley ('2006'), we will refer to this as the PRH 

(Podsiadlowski-Rappaport-Han) formul ation. 

Another f ormul ation (e.g. Ilben & Livid 119931 ; 
lYungelson et alJ 11994") takes the binding energy as being 
between the envelope mass and the combined mass of the core 
of the primary and the secondary star 



p - 

Cgr — 



G(Mi.e + M2)Mi.e 

2ai 



(7) 



4 
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(IJavis et al. 2008) 




Y\^. (RappaporL et al 1983) 




l^orb (observed) 

■ SDSS-PCEBs 

■ all WDMS-PCEBs 



log(Porb[days]y 



15 



10 - 



P^j. (Davis et al. 2008) 
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P(,^ (RappaporL et al 1983) 





J-^orb (observed) 

■ SDSS-PCEBs 

■ all WDMS-PCEBs 



il.l, I , I , ■, , 

1 3 3 

Po.b[days] 



w ■ 



Figure 1. Observed, present-day (bottom) and reconstructed post-CE (middle and top) orbital period distributions (left: logarithmic 
scale, right: linear scale) for two different prescriptions of magnetic braking. The well-defined sample of SDSS-PCEBs is shown 
in black while the gray distribution represents the entire WDMS PCEB population (IK Peg is not present in the right panel of this 
figure due to its long period compared with the rest of the sample). There is no significant difference between the two populations. 
The observed distributions as well as the reconstructed zero-age PCEB distributions show a strong peak at ~ 8 hrs and a secondary 
peak at ~ 17 hrs. The reconstructed distributions for both prescriptions of disrupted magnetic braking are very similar because most 
PCEBs are relatively young and/or contain fully convective secondary stars. 



and the initial orbital energy as the orbital energy between the 
core of the primary and the secondary at the initial binary sepa- 
ration 



■^orbj 



1 GMi,eM2 

2 a\ 



(8) 



We will refer to this as the ILY (Iben-Livio-Yungelson) formu- 
lation. 

Finally there is another scheme, use d in the binary star e volu- 
tion (hereafter BSE) code presented bv lHurlev et al.l (l2002h . that 
takes the gravitational energy in the same way as in the PRH for- 
mulation (Eq.|5]) and the initial orbital energy as in the ILY for- 
mulation (Eq.[8]l. We will refer to this as the BSE formulation. 
We compare the results obtained with the three formulations in 
Sect.lH 



5. The reconstruction algorithm 

As in NT05, we determined the possible masses and radii of 
the progenitors of the WDs in all the PCEBs listed in Table lA.ll 
and IA.2I from fits to detailed stellar-evolution models. We as- 
sumed that the observed WD mass (Mwd) is equal to the core 



mass of the giant progenitor (Mi c) at the ons e t of m ass trans- 
fer and used the equations from iHurlev et al.l ^2000) to calcu- 
late the luminosities Lg an radii Rg of all giant stars that have 
exactly such a core mass. We did this for initial masses Mi of 
1.0, 1.01, 1.02, ...Mq up to the mass for which the initial core 
mass, i.e. the core mass at the end of the main sequence, is larger 
than the observed WD mass. We also included possible pro- 
genitors in the hertzprung gap (HG) with initial masses greater 
than 1.2 MQ(to ensure a convective envelope). Because we used 
equations from lHurlev et al.l ( l200C ft for the different evolutionary 
stages instead of running the code, we set mass dependent lumi- 
nosity limits for the progenitors in different evolutionary phases. 
For stars in the HG we required the luminosity to be between the 
luminosity at the top of the main sequence and the luminosity at 
the base of the first giant branch (FGB) (i.e. Ltms < Lg < Lbgb)- 
For the FGB, the luminosity should be between the luminos- 
ity at the base and at the end of the FGB phase respectively 
(i.e. Lbgb < Lg < Luei). For the early asymptotic giant branch 
(EAGB), we required the luminosity to be between the lumi- 
nosity at the base of the AGB and the luminosity of the second 
dredge-up (i.e. Lbagb < Lg < Ldu)- Finally, for the second AGB 
(SAGB, i.e. after the second dredge-up) we required the lumi- 
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nosity to be lower than the pe ak luminosity of the first thermal 
pulse according to Eq. (29) in 'i zzard et alJ (l2004l) . For all pos- 
sible progenitors we also required q - M1/M2 to be greater 
th an a critical value (qcm), neccessary to have a CE according 
to lTout et al.l (Il997h (their Eq. 33) or greater than 3.2 for a pro- 
genitor in the HG. 

After obtaining the mass and radius of a possible WD pro- 
genitor with a core mass equal to the measured WD mass, we 
assumed that the giant radius was equal to the Roche-lobe radius 
at the onset of mass transfer. Because the secondary mass is as- 
sumed to remain constant during the CE phase, this allows us to 
determine the orbital separation just before the CE phase. The re- 
maining quantities in the CE equation are then the CE efliciency 
a and the binding energy parameter A, and we can derive aA 
for each possible progenitor. In other words, from Roche geom- 
etry and the energy equation, we get one value for aA for each 
parameter set consisting of the progenitor mass, core mass {- 
current WD mass), secondary mass (= current secondary mass), 
and final orbital period {- Pce)- In this way we obtain a range 
of values for aA for each system that corresponds to a range of 
possible progenitor masses Mi . 



6. Comparing CE prescriptions 

Before discussing below what we might be able to learn from 
reconstructing the new and much larger sample of PCEBs, we 
compare here the results obtained with the three formulations 
used to describe the CE evolution (Sect.©. Each horizontal line 
in Fig.|2]represents possible values of aA for different masses of 
the progenitor for a given WD and secondary mass. As in NT05, 
the different lines for each object represent different values of 
the WD mass within O.OSMq from the best-fit value (or within 
the error given in Table lA. 11 and lA.2l in case it exceeds O.OSMq). 
Values obtained for FGB progenitors are shown in black, while 
AGB progenitors are in blue (see the electronic version of the pa- 
per for a color version). We did not find any possible progenitor 
on the HG phase. 

Solutions for most of the non-SDSS PCEBs are also given in 
NT05 (their Fig. 6). Comparing their results with those shown in 
the left panel of Fig.|2]one finds that the two reconstruction algo- 
rithms give very similar results with the only notable difference 
that we generally find slightly more solutions for a given system. 
This is because our grid of progenitor masses is finer by a factor 
of ten (step size O.OIMq instead of O.IM0). 

Comparing the three panels of Fig.|2]it becomes obvious that 
the results obtained with the PRH and BSE algorithm are almost 
identical, with aA being a little but not significantly lower for 
the BSE scheme. There are, however, significant diff'erences be- 
tween those two formulations and the ILY scheme, which gives 
by far the lowest values. This is easy to understand as the ILY 
formulation predicts much lower values for the gravitational en- 
ergy than the PRH prescription. We also note that the ILY ver- 
sion of Egr does not contain the structural parameter A. Hence, 
we are in fact plotting a for this formulation. In general, it is dif- 
ficult - if not impossible - to judge which of the three algorithms 
for the initial conditions of CE evolution should be used. In any 
case, much of the physics is contained in the parameters a and A. 
As most calculations presented in the literature are based on the 
PRH or the BSE formalism, we will use the BSE formulation in 
the following sections to facilitate the comparison of our results 
with those obtained by other authors. 



7. The binding energy of the envelope 

The structural parameter A has generally been taken as a constant 
(typically ~ 0.5). Detailed stellar models taking into account the 
structure of the envelope show that this is a reasonable assump- 
tion as long as the internal energy of the envelope is ignored. In 
this case one o b tains A ~ 0.2 - 0.8. Howev e r, acco rding to e.g. 
iDewi & TaurisI (l2000l) : iPodsiadlowski et all (l2003h . A = 0.5 is 
not a very realistic assumption if a fraction of the internal energy 
of the envelope supports the process of envelope ejection. In this 
case, especially the extended envelopes of luminous AGB stars 
can be very loosely bound, i.e. reaching values of At 10. This is 
mainly due to the recombination-energy term. It is still not en- 
tirely clear whether this energy indeed contributes to unbind the 
envelope of the donor or if it is entirely radiated away (see e.g. 
ISoker & Harpazll2003t iHan et al. ] |2003[ for further discussion). 
However, the internal energy of the envelope might be a very im- 
portant factor to explain the e xistence of long orbital period sys- 
tems, and we therefore follow lHan et al ] (11995), who introduced 
a parameter ffth (between and 1) to characterize the fraction 
of the internal energy that is used to expell the CE. Using this, 
and calling the parameter aim (as it includes not only the thermal 
energy, but also the radiation and the recombination energy), the 
equation for the standard a-formalism becomes 

QTorbA^orb - -Bgr - aintt/int- (9) 

Alternatively one can revise A to incorporate the internal en- 
ergy t/int. If a fraction of the internal energy contributes to ex- 
pelling the envelope, the binding energy writes 

r'^' I Gm \ 
-Ebind = — — +aiMUM{m)]dtn. (10) 
Jm,, \ i-Cn) j 

Detailed calculations of this expression have been performed by 
vario us authors (e.g. Dewi & Tauris 2000; Podsiadlowski et aTj 
I2003h who demonstrate that the binding energy depends sig- 
nificantly on the mass of the giant, its evolutionary state, and 
of course, o-int. Clearly, to include the effect of the internal en- 
ergy and the structure of the envelope in the simple energy 
equation (Eq.O one may equate Ebmd with the parametrized 
binding energy from Eq. Q, keeping A variable. The latest 
version of the BSE code includes an algorithm that computes 
A in this wa y. -Ebind has been calculated using detailed stellar 
models from iPols et all (1199 8') and approximated with analyt- 
ical fits (Pols, private communication). Using this algorithm A 
is no longer a constant but depends on the mass, the evolution- 
ary state of the mass donor, and on the fraction of the internal 
energy used to expel the envelope, i.e. aint. Note that the exact 
definition of the core radius that separates the ejected envelope 
from the condensed core region in the primary is of major im- 
porta n ce for high-mass progenit ors on the FGB (iTauris & Dewil 
Ivan derSluvs etal]|2006h . As we mostly find low-mass 
progenitors on the FGB the exact definition of the core radius 
(and hence of the core mass) can be assumed to be of minor im- 
portance here. The prescription of A used in this work is based 
on the core-envelope boundary being defined as the mass shell 
where the hydrogen mass fraction becomes less than 10%. 

In the next sections we assume the efficiency of using the 
internal energy of the envelope and the orbital energy to expell 
the envelope to be equal, i.e. we use values of A that include a 
fraction o-int = aoit = Q'- Hence, the given values of a should 
be interpreted as the fraction of the total energy that is used to 
expell the envelope, independent of whether this energy has to be 
transferred from the orbit to the envelope or was already present 
in the envelope as internal energy. 
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Figure 2. Reconstructed values of a A for the three different versions of the energy equation. Black lines are progenitors in the FGB, 
and blue are for progenitors in the AGE. The results obtained with the PRH and BSE formulations (left and center) are almost 
identical and significantly higher than those obtained with ILY formulation (right panel) because in the latter case the binding 
energy at the onset of CE evolution is assumed to be significantly smaller. 



In Fig.|3]we plot the possible values of a for each PCEB in 
our sample assuming A - 0.5 (left), calculating A with the BSE 
algorithm but without internal energy (center), and including a 
fraction ai^i - a of the internal energy (right). Our results indi- 
cate that the structural parameter is quite well approximated by 
assuming A - 0.5 for FGB progenitors (in black) because there 
is hardly any difference between the black lines in the three pan- 
els. However, the effect of calculating A and including the in- 
ternal energy is of utmost importance for AGB progenitors: the 



blue lines move towards lower values of a especially if a frac- 
tion of the internal energy is assumed to contribute to the energy 
budget of CE evolution. The effect is most obvious for IK Peg 
because we only find a solution with < a < 1 if the internal 
energy is included. This result perfectly agrees with iDavis et al.l 
(.2010.) . In addition, the internal energy becomes important espe- 
cially for long orbital period systems - exactly as suggested by 
IWebbinkI (120081) . 
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Figure 3. Reconstructed values of a for all the possible progenitors of the PCEBs in our sample with A - 0.5 (left), A calculated 
using the BSE code without internal energy (center), and with A calculated including a fraction aint = o- of the internal energy 
(right). Colors are the same as in Fig.|2l While A - 0.5 seems to be a reasonable assumption for most of the FGB progenitors, 
calculating A and particularly including internal energy becomes important for progenitors on the AGB (blue). While a is only 
slightly moved towards lower values in the central panel, taking into account the internal energy leads to dramatically lower values 
of a for AGB progenitors (right panel). For example, we only find solutions for IK Peg in the range < a < 1 if a fraction of the 
internal energy is assumed to contribute to the energy budget. The vertical lines in the right panel correspond to a = 0.2 and 0.3. 



Inspecting the right panel of Fig.[3]in more detail, it becomes 
obvious that including the internal energy allows us to find so- 
lutions for all the systems in a small range of CE efficiencies, 
i.e. a - 0.2 - 0.3 (vertical lines). The upper limit of this range 
(a = 0.3) is defined by systems with massive WD (so they have 
progenitors in the AGB) and short orbital periods after the CE 



phase. In contrast, the lower limit (a - 0.2) is given by systems 
with FGB progenitors (i.e. those with low-mass WDs). 
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8. a versus y 

As mentioned above, NT05 used a similar algorithm to recon- 
struct the CE phase of double WDs and PCEBs. The problem 
they encountered can be summarized as follows: during the first 
CE phase of virtually all double WDs and for three alleged 
PCEBs (AYCet, SI 040 and IK Peg) the observed binary separa- 
tion is too large, requiring a physically unrealisticly high or even 
a negative efficiency. NT05 therefore proposed to use the angu- 
lar momentum conservation instead of the energy conservation 
because they find the angular momentum relation in agreement 
with the observed binary separations of double WDs and all the 
PCEBs in their sample. The alternative angular momentum algo- 
rithm for CE evolution (the so called y-algorithm) is described 
by 



M 



I.e 



Ml + M2 ' 



(11) 



where ^ is the relative change in angular momentum and 
is the relative change in mass. At first glance, the fact that all the 
unexplained double WDs and the three alleged critical PCEBs 
have reasonable solutions for y appears to be very attractive. 
Moreover, the obtained values of y cluster in a rather small range 
of values, i.e. 1.5 < y < 1.75, raising hope for a new and uni- 
versal prescription of CE evolution. However, this turned out to 
be an illusion as iWebbinkI (2008) recently showed that energy 
conservation is much more constraining the outcome of CE evo- 
lution. Indeed, a final energy state lower than the initial one re- 
quires the loss of angular momentum while the opposite is not 
necessarily true. In addition, Webbink (2008) showed that the 
ratio of final to initial orbital separation is extremely sensitive to 
y in the range of values proposed by NT05. 

Our large and representative sample of PCEBs now allows 
us to test both algorithms and to evaluate their predictive power 
In Fig.|4]we show the values of a (left) and y (right) for all possi- 
ble progenitors of the PCEBs in our sample. The binding energy 
parameter A has been calculated with the BSE code including 
internal energy. All the PCEBs in our sample can be simulta- 
neously reproduced with a in the range 0.2 - 0.3 (vertical lines 
in the left panel). The right panel shows that indeed literally all 
systems can be reconstructed with y = 1 .5 - 1 .75 (vertical lines). 

In Fig.|5]we investigate the effect of constraining a on the 
possible range of values for y and vice versa. In the left panel 
we show the values of a requesting 1.5 < y < 1.75, while on 
the right hand side we show the values of y if 0.2 < a < 0.3. 
Apparently, requesting y to lie in a small range of values does 
not very much constrain the values obtained for a. We still find 
the solutions for a covering basically the entire parameter space, 
i.e. < ff < 1. This confirms the suggestion of Webbink (2008) 
that virtually all possible configurations can be explained with 
similar values of y, which questions the predictive power of the 
new algorithm. 

In contrast to this, fixing a provides strong constraints on y. 
The values we obtain for y seem to have a clear dependency on 
the evolutionary stage of the WD progenitor. It is almost constant 
for progenitors in the same evolutionary stage, being higher for 
FGB progenitors and smaller for AGB progenitors. This finding 
has a straightforward physical interpretation: the envelope of a 
giant star is more tightly bound on the FGB and less bound on 
the AGB, where it is more expanded (especially on the second 
AGB). The value of y represents the ratio of the relative amount 
of angular momentum loss to the relative amount of mass loss. 
Hence, the different values of y may just reflect the simple fact 



that expelling a tightly (loosely) bound envelope requires to ex- 
tract more (less) angular momentum per unit mass. 

Once more, the findin gs descri b ed abo ve perfectly agree with 
the results obtained by IWebbinkI (l2008l) . i.e. we need to con- 
strain a to predict the outcome of CE evolution. In addition, 
the internal energy of the envelope seems to play an important 
role. Taking this into account in the energy equation leads to two 
classes of solutions in the angular momentum equation. 

9. Should a be constant? 

In most binary p opula tion synthesis models of WDMS (e.g. 
Willems & Kolbl 12004') bu t also of soft X-ray transients 
(Yungelson & Lasota 2008: iKiel etalj|2008h or extreme hori- 
zontal branch stars (Han e t alJ 2002h . the CE efficiency is as- 
sumed to be constant. Analyzing our sample of PCEBs consist- 
ing of WDs and low-mass main-sequence stars we find that we 
can reconstruct the evolutionary history of all systems assuming 
a constant value a ~ 0.2 - 0.3. 

An important question is now whether we should expect a 
to be constant for all types of PCEBs . First s teps exploring this 
have been made bv lPolitano & Welled (l2007h ; iDavis et alj (l2008l 
;201 0) who recently speculated that instead of being constant, a 
may depend on the mass of the secondary star or on the final 
orbital separation as spiraling-in deeper into the envelope may 
significantly affect t he efficiency of the ejection process. 

We here foUow iDavis et alJ (l2008h and evaluate the forma- 
tion probability for each possible progenitor of each PCEB in 
our sample. The number of primaries with masses in the range 
dMi is given by dN oc /(Mi )dMi where /(Mi ) is given by the 
initial mass function (IMF): 



/(Ml) 



-1.3 





0.29056M. 
0.15571Mj-2-2 

0.15571Mr" 



Ml/Mo < 0.1, 
0.1 < M1/M0 < 0.5, 
0.5 < Ml /Mo < 1.0, 
1.0 < Mi/Mq, 



(12) 



dKroupa et alj|1993l) . The probability that a binary forms with a 
certain initial orbital separation a-, is determined by 



/i(fli) — 





0.078636flr 



ai/Rs < 3 OTOi/Ra > lO*", 
3 < Ui/Ro < 10^ 



(13) 



(iDavis et al.ll2008h . The formation probability for each progeni- 
tor is then given by P{Mi,ai) = f{Mi)h(a(). 

In Fig.|6]we plot the weighted mean value of a for each sys- 
tem (colored points) versus the mass of the secondary star (left) 
and the orbital period the PCEB had at the end of the CE phase 
(right). Black vertical lines represent the full range of possible 
values of a. Again, we distinguish between progenitors in differ- 
ent evolutionary stages. Red points indicate systems with pro- 
genitors on the FGB, while blue points are for progenitors on 
the AGB. Given the uncertainties in the WD masses, some sys- 
tems have possible progenitors in more than one evolutionary 
stage. For those cases, we separately computed the average for 
the different type of progenitors. Finally, dashed horizontal lines 
indicate a - 0.2 and 0.3. There seems to be no dependence of 
a on the mass of the secondary star or on the period, but a large 
scatter around a - 0.2 - 0.3. 

This finding remains if we assume alternative initial mass 
distributions. We tested for two other probability distributions 
assuming that the masses of the binary components are corre- 
lated. We used niq-i) oc q2 and n{q2) oc q^ '^'^ , where ^2 = 
M2/M1. In both cases we obtained very similar results, i.e., a 
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Figure 4. Reconstructed values for a {left) and y {right) for the possible progenitors of our PCEB sample. The structural parameter 
A has been calculated including a fraction a^i = a of the internal energy of the envelope. On the left hand side, the vertical lines 
indicate the range of values were we find simultaneous solutions for all the systems in our sample, i.e. a = 0.2 - 0.3. On the right 
panel vertical lines show the range of simultaneous solutions for j proposed by NT05. 



large scatter and no relation between a and M2 or Pqe- Although 
there seems to be no correlation between a and the mass of the 
secondary or the final period, there is a clear relation between the 
averaged mean values of a and the evolutionary state of the pro- 
genitor Systems with FGB progenitors tend to have weighted 
mean values a > 0.3, while the obtained mean efficiencies for 
systems with AGB progenitors are much smaller, i.e. aSO. 1 . This 
is easily explained if one remembers that the internal energy be- 
comes very important for progenitors on the AGB moving the 



whole range of possible values of a towards smaller values (see 
Sect.|7]i. It is essential to recall here that the given values of 
a represent the fraction of the total energy that is used to ex- 
pell the envelope. In other words, the same fraction of internal 
and orbital energy are used, i.e. aint - a. However, one could 
also point out that the orbital energy must first be transferred 
to the envelope (presumably as thermal energy), in contrast to 
the energy already present in the envelope and that this would 
give rise to a different a for the two. Indeed, the systematically 
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Figures. Left panel: reconstructed values of a for y fixed between 1.5 and 1.75. Right panel: reconstructed values of y for a fixed 
between 0.2 and 0.3. If constraining y we still find rather broad ranges of possible values for a. In contrast, if we constrain the energy 
efficiency to be a = 0.2 - 0.3, the values of y cluster in a small range of values and there is a clear dependency on the evolutionary 
stage of the progenitor of the primary, which reflects the fact that expelling tightly bound evelopes extracts more angular momentum 
per unit mass from the binary (see text for details). 



lower weighted mean values of a for AGB progenitors may re- 
flect different efficiencies for the orbital and internal energy. If 
aint is small, the required fforb will increase especially for sys- 
tems with AGB progenitor So, an alternative to a =const. might 
be fforb =const and aint =const but aim < aorb- A detailed dis- 
cussion of this alternative possibility is beyond the scope of this 
paper though. 



As a final remark we emphasize that the weighted mean val- 
ues discussed above are lacking a physical meaning. We used 
these values here only to test for possible dependencies of a 
that are missing in the energy equation, which does not seem 
to be the case. Therefore, a =const. or at least aorb =const. and 
(Tint =const., which corresponds to the assumption that the most 
important dependencies are included in the used energy equation 
remains the currently most reasonable prescription. 
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Figure 6. Weighted mean values of a versus M2 (left) and f ce (fight). Red is for systems with FGB progenitors, while blue is for 
AGB progenitors. The full range of possible values of a is given by the solid vertical lines. Dashed horizontal lines are for a = 0.2 
and 0.3). 



10. Discussion 

The results obtained in the previous sections can be summarized 
as follows: For all systems in our sample, the largest sample of 
one specific type of PCEBs that is currently available, we find 
possible progenitors assuming energy conservation if the inter- 
nal energy of the envelope is taken into account. For each indi- 
vidual system the possible solutions cover rather broad ranges 
of values for the CE efficiency a. However, there exists only a 
small range of values, i.e. a - 0.2 - 0.3 for which we find so- 
lutions for all the systems in our sample. This means that, if a 
universal value for the CE efficiency does exist, it should lie in 
this range. A plausible alternative to such a universal value for 
a is to assume that the fraction of the orbital energy exceeds the 
fraction of the internal energy that is used to expell the enve- 
lope, i.e. ffint < ffoib- In addition, we have shown that the energy 
budget constrains the outcome of CE evolution much more than 
the alternative angular momentum equation. In this section we 
discuss our results in the context of recent theoretical and ob- 
servational results in the field of close compact binary formation 
and evolution. 



10.1. Hydrodynamical simulations 



Soon after iPaczvnskil (1 1 976h outlined the basic ideas of CE evo- 
lution, the first hydrodynamical simulations in one dimension 
were carried out (Taam et al. 1978; Mever & Mever-Hofmeiste^ 
Il979h . Based on these early studies two and three dimen- 
sional models have been de velope d in the last decades 



(e.g. Bodenheimer & TaamI 11984 iTaam & BodenheimeJ 



1989"; Sandauist et al.' '2000*). For a recent review see e.g 
^Taam & Ricker (2006). The most important findings of hydro- 
dynamical simulations of CE evolution are perhaps the relatively 
short duration of CE evolution (SlOOOyrs) and the preference of 
ejecting matter in the orbital plane. In addition, as most particles 
are predicted to leave the CE with velocities exceeding the 
minimum escape speed, the predicted CE efficiency is less than 
40 - 50%, i.e. ai:Q.4 - 0.5. This result agrees quite well with our 
finding of a = 0.2 - 0.3. However, one should note that current 



hydrodynamical simulations still cannot follow the entire CE 
evolution basically because of the large ranges of timescales and 
length scales that have to be numerically resolved. Therefore, 
even the most detailed hydrodynamical simulations still have to 
be considered as rather rough approximations. 

10.2. Binary population synthesis 

An alternative way to constrain the CE efficiency is to perform 
binary population studies and compare the predictions with the 
observed properties of PCEBs. These simulations have become 
popular in last 10-20 years and have been carried out for a large 
variety of different PCEB populations. We here briefly review 
the main results. 



10.2.1. WDMS binaries 



Th e pop ulation of WDMS binaries has been first simulated by 
de Kooil dl992h and Ide Kool & Ritted dl993h . Ide Kool & Ritted 
(1 19931) incorporated observational selection effects to compare 
their predictions with the - very small and biased - observed 
populations they had at hand. Interestingly, for a - 0.3 and 
M2 randomly taken from the IMF they predict PCEB orbital 
period distributions rather similar to the observed distribu- 
tion (see Sec t.[3l). However, the selection effects applied by 
Ide Kool & Rit ter ( 1993) have been designed for blue color sur- 
veys such as the Palomar Green survey and are not applicable 
to our new SDSS PCEB sample. In addition, one should take 
into account that the app roximations to stellar evolution used by 
Ide Kool & Rifted d 19931) have been much cruder than the models 
that are available today and that they did not include the internal 
energy of the envelope. 

An update of this early work was carried out by 

IWiUems & Kolbl (|200l. using more detailed analytical fits to 
stellar evolution f Hurlev et afll2000l) . Their PCEB orbital pe- 
riod distribution peaks at about one day, i.e. at a significantly 
longer pe riod than the observed sample. However, one should 
note that IWiUems & Kolbl d2004l) computed formation models 
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for PCEBs, but did not follow the subsequent angular momen- 
tum loss by magnetic braking and gravitational radiation. In ad- 
dition, no observational biases are incorporated in their predi- 
tions. Hence, we advocat e caution when comparing the predic- 
tions of IWillems & Kolbl (|2"004.) with observed samples. 

Full b inary population studies of PC EBs have been per- 
formed b V 'Po litano & W eiler ('20061 l2007l) . They tested different 
formulations of a and discussed the influence on the predicted 
distributions. The resulting orbital period distributions peak at 
f orb ~ 3 days and the overall shape does not change signifi- 
cantly for different prescriptions of the CE efficiency. Again, as 
observational selection effects have not been incorporated, it is 
difficult, if not impossible, to compare the predicted distributions 
with the measured orbital perio d distri butions shown in Fig.[T] 

Most recently, iDavis et aP (12010 ) published a work pre- 
senting comprehensive population synthesis studies of PCEBs. 
Perhaps most importantly, for the first time the PCEB population 
has been simulated including variable values of A. Comparing 
their predictions with the observations, iDavis et al ] (120101) find 
a disagreement in the orbital period distributions, i.e. the pre- 
dicted distributions peak at f orb ~ 1 day declining smoothly at 
longer periods, while observ ations indica t e a rat her steep decline 
at Poib ~ 1 day. However, iDavis et aTl (1201 Ol) compared their 
predictions with a small sample of PCEBs identified through 
various detection channels. Thanks to our concentrated follow- 
up of WDMS binaries from SDSS, the number of known PCEBs 
has increased by more than a factor of two, and this new SDSS 
PCEB sample is less aff'ected by observational biases (Gansicke 
et al. 2010, MNRAS in p rep). In addition, the parameter space 
explored by Davis et aP (1201 Oh is still rather small. While the 
CE efficiency has in general been varied over a wide range of 
values (q' = 0.1 - 1), only one model with variable values of A 
assuming a = 1 has b een calculated. Finally, one should keep in 
mind that D avis et al. (|2010) interpolated the tables provided by 
iDewi & TaurisI (l200ol) to determine A, which probably leads to 
underestimating A for large radii. 

We conclude that binary population synthesis (BPS) simula- 
tions using a - aint = 0.2 - 0.3 and including a proper treat- 
ment of A do not yet exist. Hence, it might not be too surprising 
that predicted period distributions disagree with the observation. 
Reducing a and incorporating the internal energy should lead to 
predicting less systems with Porb ^ 1 day. Therefore we antici- 
pate that applying our results may bring theory and observations 
into agreement. In addition, the next generation of BPS simula- 
tions should take into account observational biases as detailed as 
possible. The importance of this might be indicated by the basic 
agreement between the predictions by Ide Kool & Ritted (Il993h 
and our observed sample. 

10.2.2. Extreme horizontal branch stars 

Extreme horizontal branch stars (EHB, also known as hot sub- 
dwarfs) are hehum-burning stars with v ery thin hydrogen en- 
velopes (iHeber et al.l Il986i ISaffer et al.l [l994.) . To explain the 
formation of these stars several scenarios have been discussed 
mostly based on single-star evolution (e.g. lKilkennv et ani997t 
[Green et al. 1986). However, as most EHB stars appear to be 
members of clos e binary s y stems , the b inary-formation chan- 
nel proposed by iHan et all (l2002i l2003h has become a popu- 
lar alternative. These authors favored a rather high efficiency 
(a ~ 0.75) when compared to the value we obtain from our 
sample. However, one should note that Han et al. (2003) did not 
explore the full parameter space and did not generally exclude 
lower values of a. 



An interesting option to constrain a might be to mea- 
sure the binary fraction of EHB stars in globular clusters. 
Moni Bidin et al. (2006) find the binary fraction in NGC 6752 to 
be much lowe r (~ 4%) tha n in the field (270%). A s spe c ulated 
by Moni Bidin et al.l and confirmed later bv lHanI (l2008h 

and Moni Bidin et al. (2009) this can be explained within the 
binary-formation scenario, as the binary fraction among EHB 
st ars in clusters is expected to decrease with time. According 
to iHanl (120081) the binary fraction - age relation is rather sen- 
sitive to the assumed CE efficiency. First results seem to favor 
high values of a. However, significantly more measurements of 
the binary fractions among EHB stars in globular clusters are 
required to derive clear constraints. 

10.2.3. Low-mass X-ray binaries 

The efficiency of CE evolution is of outstanding importance in 
the context of compact binaries descending from more massive 
stars too. For example, the existence of low-mass X-ray binaries 
(LMXBs) in our galaxy has been difficult to explain within the 
CE picture as low-mass companions appea r to be unable to un- 
bind t he envelope of a massive primary star dPodsiadlowski et aH 
'2003^ and one therefore expects most systems to merge instea d 
of forming a LMXB. As shown by Podsiadlow ski et all (l2003h . 
the predicted formation rate of LMXBs is much lower than in- 
dicated by observations even for a - 1. This is explained by 
the huge binding energy of envelopes ar ound massive cores, i . 
e. A<0.1. As a solution for this problem. I xiel & HurlevI (l2006h 
proposed a reduced mass-loss for helium stars and brought into 
agreement binary pop ulations synthesis and o bservations for 
a ~ 1.0 (but see also lYungelson & Lasotal 12008) . In any case, 
current models seem to be unable to reproduce the observed pop- 
ulation of LMXBs assuming a rather low value of or = 0.2 - 0.3 
as we find for our sample of PCEBs. This indicates that either 
the efficiency is different for LMXBs or that the uncertainties 
in evolutionary models of very massive late AGB stars strongly 
affect the predictions of BPS. 



11. Conclusion 

We have developed a new algorithm to reconstruct CE evolution 
of PCEB stars. We included a proper treatment of the binding 
energy parameter A taking into account the internal energy of 
the envelope. We have applied the new algorithm to the largest 
and most homogeneous sample of PCEBs currently available. 
The basic result of this investigation can be summarized with 
the following four statements: 

- A reasonable prescription of the CE evolution of PCEBs con- 
taining a WD primary and late-M spectral-type secondary is 
given by the energy equation if the internal energy of the en- 
velope is included. 

- The energy equation is much more constraining the outcome 
of CE evolution and the predictive power of the angular mo- 
mentum equation is limited. 

- If there is a universal value of a, it must be in the range of 
0.2-0.3. 

- There are no indications for a dependence of a on the mass 
of the secondary star or the orbital period. 

Despite these findings, it is still unclear whether a univer- 
sal constant value of a can explain CE evolution in general. 
Answering this question requires to observationally establish 
representative and large samples of all types of PCEBs, i.e. not 
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only WDMS, but also neutron star/black hole PCEBs. However, 
if such a value exists, our result of a = 0.2 - 0.3 can be inter- 
preted as a definitive answer to one of the important questions in 
close compact binary evolution, especially as there seems to be 
no dependence of the CE efficiency on the mass of the secondary 
star or the final orbital period. 
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Table A.l. Properties of the SDSS PCEBs 



Object 






P 
(d) 


MwD 

(Mo) 






M2 

(Mo) 




(K) 


(Gyr) 


^CE 

(d) 


Ref. 


SDSS 1435 


+ 


3733 


0, 


.126 


0.505 


+ 0, 


.025 


0. 


.218 


+ 


0. 


.028 


12392 


0.275 


0.133 


1 


SDSS 1648 


+ 


2811 


0, 


.131 


0.630 


+ 0, 


.520 


0. 


.320 


+ 


0, 


.060 


13432 


0.284 


0.142 


2 


SDSS0052 


- 


0053 


0, 


.114 


1.220 


+ 0, 


.370 


0. 


.320 


+ 


0. 


.060 


16111 


0.421 


0.143 


3 


SDSS2123 


+ 


0024 


0, 


.149 


0.310 


±0, 


.100 


0. 


.200 


+ 


0. 


.080 


13279 


0.000 


0.149 


4 


SDSS1529 


+ 


0020 


0, 


.165 


0.400 


±0, 


.040 


0. 


.260 


+ 


0. 


.040 


14148 


0.300 


0.170 


3 


SDSS1411 


+ 


1028 


0, 


.167 


0.520 


±0, 


.110 


0. 


.380 


+ 


0. 


.070 


30419 


0.009 


0.188 


4 


SDSS 1548 


+ 


4057 


0, 


.185 


0.646 


±0, 


.032 


0. 


.174 


+ 


0. 


.027 


11835 


0.416 


0.191 


1 


SDSS0303 


- 


0054 


0, 


.134 


0.912 


+ 0, 


.034 


0. 


.253 


+ 


0. 


.029 


;8000 


;2.24 


;0.20 


1 


SDSS2216 


+ 


0102 


0, 


.210 


0.400 


+ 0, 


.060 


0. 


.200 


+ 


0. 


.080 


14200 


0.297 


0.212 


4 


SDSS 1348 


+ 


1834 


0, 


.249 


0.590 


±0, 


.040 


0. 


.319 


+ 


0. 


.060 


15071 


0.184 


0.251 


5 


SDSS0238 


- 


0005 


0, 


.212 


0.590 


±0, 


.220 


0. 


.380 


+ 


0. 


.070 


21535 


0.045 


0.261 


4 


SDSS2240 


- 


0935 


0, 


.261 


0.410 


±0, 


.080 


0. 


.250 


+ 


0. 


.120 


12536 


0.443 


0.263 


4 


SDSS 1724 


+ 


5620 


0, 


.333 


0.420 


+ 0, 


.010 


0. 


.360 


+ 


0. 


.070 


35746 


0.000 


0.333 


3 


SDSS2132 


+ 


0031 


0, 


.222 


0.380 


±0, 


.040 


0. 


.320 


± 


0. 


.010 


16336 


0.179 


0.333 


4 


SDSSOllO 


+ 


1326 


0, 


.333 


0.470 


±0, 


.020 


0. 


.310 


+ 


0. 


.050 


25167 


0.051 


0.333 


1 


SDSS1212 


- 


0123 


0, 


.333 


0.470 


+ 0, 


.010 


0. 


.280 


+ 


0. 


.020 


17304 


0.191 


0.334 


6 


SDSS1731 


+ 


6233 


0, 


.268 


0.450 


±0, 


.080 


0. 


.320 


± 


0, 


.010 


16149 


0.228 


0.361 


4 


SDSS 1047 


+ 


0523 


0, 


.382 


0.380 


+ 0, 


.200 


0. 


.260 


+ 


0, 


.040 


12392 


0.417 


0.384 


7 


SDSS 1143 


+ 


0009 


0, 


.386 


0.620 


+ 0, 


.070 


0. 


.320 


+ 


0. 


.010 


16910 


0.138 


0.411 


4 


SDSS2114 


- 


0103 


0, 


.411 


0.710 


±0, 


.100 


0. 


.380 


+ 


0. 


.070 


28064 


0.018 


0.416 


4 


SDSS2120 


- 


0058 


0, 


.449 


0.610 


+ 0, 


.060 


0. 


.320 


+ 


0. 


.010 


16149 


0.156 


0.450 


4 


SDSS 1429 


+ 


5759 


0, 


.545 


1.040 


+ 0, 


.170 


0. 


.380 


± 


0, 


.060 


16336 


0.401 


0.566 


5 


SDSS 1524 


+ 


5040 


0, 


.590 


0.710 


+ 0, 


.070 


0. 


.380 


+ 


0. 


.060 


19640 


0.109 


0.601 


5 


SDSS2339 


- 


0020 


0, 


.655 


0.840 


+ 0, 


.360 


0. 


.320 


+ 


0. 


.060 


13266 


0.508 


0.657 


3 


SDSS1558 


+ 


2642 


0, 


.662 


1.070 


+ 0, 


.260 


0. 


.319 


+ 


0. 


.060 


14560 


0.609 


0.664 


5 


SDSS1718 


+ 


6101 


0, 


.673 


0.520 


+ 0, 


.090 


0. 


.320 


+ 


0. 


.010 


18120 


0.075 


0.678 


4 


SDSS1414 


- 


0132 


0, 


.728 


0.730 


±0, 


.200 


0. 


.260 


± 


0, 


.040 


13904 


0.329 


0.729 


7 




4- 


nn/1 1 


0, 


.728 


0.900 


±0, 


.1 jU 


u. 


.380 


+ 


0. 


,U1U 


lOJ /Z 




U. 1 jy 


J 


SDSS 1705 


+ 


2109 


0, 


.815 


0.520 


+ 0, 


.050 


0. 


.250 


+ 


0. 


.120 


23613 


0.023 


0.815 


4 


SDSS 1506 




0120 


1, 


.051 


0.430 


±0, 


.130 


0. 


.320 


± 


0, 


.010 


15422 


0.251 


1.057 


4 


SDSS1519 


+ 


3536 


1, 


.567 


0.560 


±0, 


.040 


0. 


.200 


+ 


0. 


.080 


19416 


0.065 


1.567 


4 


SDSS 1646 


+ 


4223 


1, 


.595 


0.550 


+ 0, 


.090 


0. 


.250 


+ 


0. 


.120 


17707 


0.093 


1.595 


4 


SDSS0924 


+ 


0024 


2, 


.404 


0.520 


+ 0, 


.050 


0. 


.320 


+ 


0. 


.010 


19193 


0.059 


2.404 


4 


SDSS2318 




0935 


2, 


.534 


0.490 


+ 0, 


.060 


0. 


.380 


± 


0. 


.070 


22550 


0.026 


2.534 


4 


SDSS 1434 


+ 


5335 


4, 


.357 


0.490 


+ 0, 


.030 


0. 


.320 


± 


0, 


.010 


21785 


0.030 


4.357 


4 



D lPvrzaset al. (2009); (2) Pyrzas et al. (2010, in prep.); (3) Rebassa-Manseraas et al. (2008); (4) Giinsicke et al. (201 0, MNRAS in 
prep.); (5) Nebot-Gomez-Moran et al. (2010, in prep.); iNebot Gomez-Moran et al.i t2009.) ; (7) iSchreiber et alj J2003) . 
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Object 



Alt. Name 



P 
(d) 



MwD 

(Mo) 



Mo 
(Mo) 



(K) 



(Gyr) 



(d) 



Ref. 



WD0137-3457 

GD448 

NNSer 

LTT560 

MS Peg 

LM Com 

CC Cet 

CSS 080502 

RRCae 

BPM6502 

GKVir 

EC 14329 - 1625 
EC 12477 - 1738 
DECVn 

EC 13349 - 3237 

RXJ2130.6 + 4710 

HZ 9 

UXCVn 

UZSex 

EGUMa 

REJ2013 +4002 

REJ1016-0520 

INCMa 

Feige 24 

IK Peg 



HRCam 
PG1550+ 131 

GD245 

PG 1224 + 309 

PG 0308 + 096 

LFT 349 
LTT3943 
PG 1413 + 015 



J 1326 + 4532 



HZ 22 

PG 1026 + 0014 
Case 1 



J0720-3146 

FS Cet 

BD + 18 4794 



0.080 
0.103 
0.130 
0.148 
0.174 
0.259 
0.284 
0.149 
0.303 
0.337 
0.344 
0.350 
0.362 
0.364 
0.470 
0.521 
0.564 
0.570 
0.597 
0.668 
0.706 
0.789 
1.260 
4.232 
21.722 



0.390 
0.410 
0.535 
0.520 
0.480 
0.450 
0.390 
0.350 
0.440 
0.500 
0.510 
0.620 
0.610 
0.530 
0.460 
0.554 
0.510 
0.390 
0.650 
0.640 
0.560 
0.600 
0.570 
0.570 
1.190 



+ 0.035 
± 0.010 
+ 0.012 
± 0.120 
± 0.020 
± 0.050 
± 0.100 
± 0.040 
+ 0.022 
± 0.050 
± 0.040 
± 0.110 
± 0.080 
+ 0.040 
+ 0.110 
+ 0.017 
+ 0.100 
± 0.050 
± 0.230 
± 0.030 
± 0.030 
± 0.020 
± 0.030 
± 0.030 
± 0.050 



0.053 : 
0.096 : 
0.111 : 

0.190: 

0.220 : 
0.280 : 

0.180: 
0.320 : 
0.182: 
0.170: 
0.100: 
0.380 : 
0.380 : 
0.410 : 
0.500 : 

0.555 : 

0.280 : 
0.420 : 
0.220 : 
0.420 : 
0.180: 
0.150: 
0.430 : 
0.390 : 
1.700: 



0.006 
0.040 

: 0.004 

0.050 

: 0.020 
: 0.050 
: 0.050 
: 0.000 

0.013 
0.010 
0.000 
0.070 
0.070 
0.060 
0.050 
0.023 
0.040 
0.000 
0.050 
0.040 
: 0.040 
: 0.020 
: 0.030 
: 0.020 
0.100 



16500 
19000 
57000 
7500 
22170 
29300 
26200 
17505 
7540 
21000 
48800 
14575 
17718 
8000 
35010 
18000 
17400 
28000 
19900 
13100 
49000 
55000 
52400 
57000 
35500 



0.179 
0.118 
0.001 
1.040 
0.027 
0.032 
0.000 
0.130 
2.037 
0.036 
0.002 
0.220 
0.113 
0.895 
0.000 
0.088 
0.086 
0.000 
0.084 
0.313 
0.002 
0.002 
0.002 
0.001 
0.027 



0.082 
0.104 
0.130 
0.162 
0.174 
0.259 
0.284 
0.288 
0.313 
0.337 
0.344 
0.421 
0.402 
0.518 
0.470 
0.530 
0.564 
0.570 
0.597 
0.688 
0.706 
0.789 
1.260 
4.232 
21.722 
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4 
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